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Abstract. We present spectroscopic observations of a large sample of Herbig Ae stars in the 10/rm spectral 
region. We perform compositional fits of the spectra based on properties of homogeneous as well as inhomogeneous 
spherical particles, and derive the mineralogy and typical grain sizes of the dust responsible for the 10 pm emission. 
Several trends are reported that can constrain theoretical models of dust processing in these systems: i) none of 
the sources consists of fully pristine dust comparable to that found in the interstellar medium, it) all sources 
with a high fraction of crystalline silicates are dominated by large grains. Hi) the disks around more massive 
stars (M>2.5M0, L>6OL0) have a higher fraction of crystalline silicates than those around lower mass stars, 
tv) in the subset of lower mass stars (M<2.5M0) there is no correlation between stellar parameters and the 
derived crys tallinity of the dust. Th e correlation between the shape and strength of the 10 micron silicate feature 
reported bv ivan Boekel et alJ J2nn3h is reconfirmed with this larger sample. The evidence presented in this paper 
is combined with that of other studies to present a likely scenario of dust processing in Herbig Ae systems. We 
conclude that the present data favour a scenario in which the crystalline silicates are produced in the innermost 
regions of the disk, close to the star, and transported outward to the regions where they can be detected by means 
of 10 micron spectroscopy. Additionally, we conclude that the final crystallinity of these disks is reached very soon 
after active accretion has stopped. 


Key words. Stars: circumstellar matter: pre-main- 
sequence stars - Infrared: ISM: lines and bands 

1. Introduction 

A characteristic shared by many young, low and interme¬ 
diate mass stars is the presence of a strong infrared excess. 
This radiation emerges from circumstellar dust grains left 
over from the star formation process. The material is be¬ 
lieved to reside in a disk which is formed as a result of 
angular momentum conservation in the collapsing molec¬ 
ular cloud. After an initial strong accretion phase, a much 
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longer pre-main-sequence phase ensues during which the 
disk slowly dissipates, and possibly planets are formed. 

The interstellar dust which finds its way into a proto¬ 
planetary disk will undergo large changes in average size 
and chemical composition. These changes trace the pro¬ 
cess of disk dissipation and planet formation. Our solar 
system contains a precious record of the processes that 
took place during its formation. Comparison of this record 
to what is observed around young pre-main-sequence stars 
provides important insight into the history of our own so¬ 
lar system, and helps to constrain planet formation mod¬ 
els. 

The infrared spectral region is rich in vibrational reso¬ 
nances of abundant dust species. Therefore, infrared spec¬ 
troscopy can be used to determine the composition of 
dust in proto-planetary disks, as well as constrain the 
size and shape of the dust grains. Analysis of the in¬ 
frared (IR) dust emission features originating from the 
disk surface layer can be used to establish to what ex¬ 
tent the dust composition in the disk has evolved away 
from that seen in the interstellar medium (ISM). For 
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insta nce, crystalline silicates appear absent in the ISM 
(e.g. iDemvk et alJ l200nl: iKemper et al lill but are a 
substantial component i n (some) comets and in inter¬ 
planetary dust part icles llMacKinnon fc Rietmeiiedll987t 
iBradlev et alJll992h . Clearly, the refractory material in 
the proto-solar cloud went through large changes as the 
solar system was formed. It should be kept in mind that 
spectroscopy in the 10 micron region is sensitive to a lim¬ 
ited grain size range: large grains (with sizes above a few 
/rm, depending on chemical composition and wavelength) 
show only weak spectral structure and do not contribute 
significantly to the infrared emission features. 

In this work we study the composition of dus t in the 
circu mstellar environment of Herbig Ae/Be stars llHerbid 
EmS), using infrared spectroscopy. We restrict our study 
to a sub-group of mostly late B and A-F type stars (here¬ 
after HAe stars). These stars show little or no optical ex¬ 
tinction and low mass accretion r ates, as derived from ra¬ 
dio analysis jSkinner et al.lll993jl. and the lack of signifi¬ 
cant veiling in optical spectra. For these lower mass mem¬ 
bers of the Herbi g class, evidence for the disk hypothesis 
is compelling (e.g .lMannings fc SargentllQQ' ^lGradv et alJ 
l200l[ lAiigereau et alJ 200l( lEisner et al.ll20(l^ 

The observed SEDs of these HAe stars can very well 
be explained with models for gas-rich, passively heated 
disks in hydrostatic equ ilibrium and a puffed-up inner rim 
toullemond et ani200lll . These models indicate that the 
emission observed at near-IR wavelengths is dominated by 
the inner rim, while the mid-IR spectrum has a large con¬ 
tribution from the dust grains in the warm surface layers 
of the disk, typically located between a few to several tens 
of Astronomical Units (AU) from the star. Therefore, mid- 
IR spectroscopy predominantly relates to the comp osition 
of the dust in this surface layer. Van Boekel et al. il2nn.‘ih 
argue that due to the turbulent nature of the disks, the 
small grains observed at the disk surface are well coupled 
to those in the disk mid-plane, and so the observations 
of the surface layers likely bear relevance for the overall 
small grain population. 

Observations of the dust composition in HAe disks 
have revealed a very rich mineralogy, and strong source 
to source variations of the dust composition (e.g. 
iBouwman et al ] |200ll) . Some stars show strong 9.7 and 
18 //m amorphous silicate emission, with a band shape 
and strength very similar to that seen in the interstellar 
medium. Other objects have only weak silicate emission, 
with sub-structure at 9.2, 10.6 and 11.3 ^m due to crys¬ 
talline silicates. At longer wavelengths, high-quality ob¬ 
servations only exist for a small number of stars, so far all 
obtained with the Infrared Space Observatory (ISO). This 
situation is rapidly improving due to ongoing observations 
with the Infrared Spectrograph on board of the Spitzer 
Space Telescope. The ISO spectra reveal a mineralogy 
dominated by the crystalline silicates fors terite and en- 
statite, i.e. Mg-rich, Fe-poor materials (e.g. iMalfait et all 

ll998htlBouwman et al.l20(HtlMeeus et al.l200ll hereafter 

MEOl). A small group of stars lacks silicate emission, 
but shows prominent emission from Polycyclic Aromatic 


Hydrocarbons (PAHs) at 3.3, 6.2, 7.7-7.9, 8.6 and 11.3 ^m. 
Many stars show a combination of silicates and PAHs. 

In recent years, several investigators have attempted 
to find correlations between the properties of the dust 
in the disk of Herbig stars on the one hand, and global 
properties of the disk (e.g. disk geometry) and/or the star 
(e.g. ma ss, luminosity, age, bin a rity) on the other hand 
fMEO l: iBouwman et al. l200lt lAcke fc van den Anckeil 
IS). Perhaps the most promising results so far are re¬ 
lations between the dust properties and the shape of the 
SEDs of the disk; for instance, the PAH bands are on av¬ 
erage stronger in sources with SEDs that rise at far-IR 
wavelengths (MEOl). This has been interpreted in terms 
of the disk geometry in the following way. Relatively red 
SEDs correspond to flaring disks, that have a large sur¬ 
face which is directly irradiated by the star. If PAHs 
are present in this flaring disk surface layer, they will 
contribute to the e mission in the familiar PAH bands. 
Van Boekel et al. ll2nn4hll studied the spatial distribu¬ 
tion of the PAH emission in HD 97048, which was found 
to be extended on a scale of l-2arcsec (~250AU) but 
clearly associated with the disk. It is likely that also in 
the other stars of our sample, the PAH emission origi¬ 
nates from the outer disk region. Disks that lack a flaring 
outer region will have much les s prominent PAH emission. 

I A eke fc van den Anckeil (|2nn4ll confirm this relation be¬ 
tween PAH band strength and S EP shape using a s ample 
of about 50 Herbig Ae/Be stars. IXcke et all ll2004ll note 
that stars with flat far-IR SEDs on average have flat mil¬ 
limeter spectral slopes, suggesting that the cold mid-plane 
grains in these sources have grown to larger size than the 
corresponding grains of stars with rising far-IR SEDs. 

Despite considerable efforts, it has proven difficult so 
far to determine relations between stellar and dust pa¬ 
rameters. This has prompted us to carry out a compre¬ 
hensive spectral survey at 10 /rm of isolated HAe stars. 
Our goal is to establish relations between star and dust 
properties that are of relevance for our understanding of 
the evolution of dusty disks around young stars, by in¬ 
creasing the number of stars for which mid-IR spectra 
are available. We have used t he Thermal Infrared M ulti 
Mode Instrument 2 fTIMMI2. iHeima.nn et a,1]l^99S^ . at¬ 
tached to the 3.6 m telescope of the European Southern 
Observatory for our spectral survey. All known optically 
bright HAe stars accessible from the La Silla observatory 
were included in our initial list of targets. In practice, high 
quality data could be obtained for stars with 10 ^m fluxes 
of about 3 Jy or more. 

Here we present the results of our spectroscopic survey 
at 10/im. In total, we obtained high quality spectra of a 
sample of 24 HAe stars, introduced in section 13 We re¬ 
port on the observations and data reduction in section |3 
An overview of the spectra is given in section 01 In sec¬ 
tion [3 we present compositional fits to the silicate feature 
observed in most of the sources, using the optical proper¬ 
ties of minerals commonly found in circumstellar disks. In 
sectionElwe discuss the implications of our results for our 
understanding of dust processing in HAe disks. 
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Table 1. Basic parameters of the sample stars. In columns 1 and 2 we give the index number by which a star can be 
identified in tables and figures throughout this work, and the name of the star, respectively. The classification of the 
sources according to MEOl is listed in the third column. In column 5 the distances, derived from direct (Hipparcos) 
parallax measurements or by association to a star forming region (SFR, column 4), are given. For the stars where 
the distance is determined by association with a SFR we assume an error of 30 % in the distance. The spectral type 
according to the MK classihcation, the effective temperature and the stellar luminosity are given in columns 6, 7 
and 8, respectively. Mass and age estimates for most stars, as derived by comparing their positions in the HR diagram 
to theoretical pre-main-sequence tracks, are given in columns 9 and 10. For HD 101412 we have no reliable distance 
estimate, and can therefore not determine its luminosity, mass and age. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

# 

star 

group 

SFR 

d 

Sp.Type 

log Teft 

log L 

Mass 

log(Age) 





[pc] 


[K] 

[Lol 

[Mq] 

[yr] 

1 

AB Aur 

la 

L519 

144tfr 

A0Ve±sh 

3.979 

1.67 

2.4±0.2 

6.3±0.2 

2 

UX Ori 

Ha 

Orion OB la 

340 ± 102 

A4IVe 

3.925 

1.68 

2.5±0.3 

6.3±0.4 

3 

HD 36112 

la 


2041“ 

A5IVe 

3.911 

1.35 

2.0±0.3 

6.5±0.3 

4 

HK ORI 

Ha 

Orion OB la 

340 ± 102 

A4pevar 

3.927 

0.87 

1.7±0.3 

>6.8 

5 

HD 245185 

la 

Orion OBla 

340 ± 102 

AOVe 

3.979 

1.26 

2.2±0.3 

>6.8 

6 

V380 Ori 

Ha 

Orion OBlc 

510 ±153 

Al:e 

3.965 

1.88 

2.8±0.5 

6.2±0.4 

7 

HD 37357 

Ha 


240 ± 72^ 

A2Ve 

4.021 

1.47 

2.4±0.4 

>7.0 

8 

HD 37806 

Ha 

Orion OB lb 

470 ± 141 

A2Vpe 

3.953 

2.13 

3.6±0.8 

5.9±0.4 

9 

HD 95881 

Ha 

Sco OB2-4? 

118 ±35 

A2HI/IVe 

3.954 

0.88 

1.7±0.2 

>6.5 

10 

HD 97048 

Ib 

Ced 111 

1751^® 

B9.5Ve±sh 

4.000 

1.64 

2.5±0.2 

>6.3 

11 

HD 100453 

Ib 


1111“ 

A9Ve 

3.869 

0.90 

1.7±0.2 

7.0±0.1 

12 

HD 100546 

la 

Sco OB2-4? 

1031^ 

B9Vne 

4.021 

1.51 

2.4±0.1 

>7.0 

13 

HD 101412 

Ha 



B9.5V 





14 

HD 104237 

Ha 

Cha III 

1161^ 

A4IVe±sh 

3.925 

1.54 

2.3±0.1 

6.3±0.1 

15 

HD 135344 

Ib 

Sco OB2-3 

140 ± 42 

F4Ve 

3.819 

0.91 

1.6±0.2 

6.9±0.3 

16 

HD 139614 

la 

Sco OB2-3 

140 ± 42 

A7Ve 

3.895 

0.91 

1.7±0.3 

>7.0 

17 

HD 142666 

Ha 

Sco OB2-2 

145 ± 43 

A8Ve 

3.880 

1.13 

1.8±0.3 

6.8±0.4 

18 

HD 142527 

la 


198l®° 

F7HIe 

3.796 

1.46 

2.5±0.3 

6.0±0.4 

19 

HD 144432 

Ha 

Sco OB2-2 

145 ± 43 

A9IVev 

3.866 

1.01 

1.8±0.2 

7.0±0.3 

20 

HD 144668 

Ha 

Lupns 3 

207131 

A5-7HI/IVe 

3.899 

1.94 

3.2±0.5 

5.7±0.3 

21 

HD 150193 

Ha 

Sco OB2-2 

1501“ 

A2IVe 

3.953 

1.38 

2.3±0.2 

>6.3 

22 

HD 163296 

Ha 


1221“ 

A3Ve 

3.941 

1.38 

2.0±0.2 

6.7±0.4 

23 

HD 169142 

Ib 

Sco OB2-1 

145 ± 43 

A5Ve 

3.914 

1.16 

2.0±0.3 

6.9±0.3 

24 

HD 179218 

la 

(L693) 

243lg 

B9e 

4.021 

2.00 

2.9±0.5 

6.1±0.4 


^ We adopted a distance of 240 pc so that the luminosity of this star matches its spectral type. 


2. The sample stars 


Our sample of stars was selected from a larger list of 
candid ate) Herbig Ae/Be stars studied bv iMalfait et al.l 
1998all . This list was constructed by comparing the 
positions of stars in the Smithsonian Astrophysical 
Observatory Star Catalog with positions in the Infrared 
Astronomical Satellite (IRAS) point source catalogue, and 
subsequent follow-up studies to find the Herbig star can¬ 
didates. Clearly, such a way to select stars may introduce 
biases. Our sample does not contain stars that are heav¬ 
ily obscured due to e.g. on-going accretion or an edge-on 
dusty disk. Instead, our sample is dominated by stars with 
a relatively “clean” environment, low optical extinction, 
and disk orientations that are not close to edge-on. An 


overview of the stars in the sample and their basic param¬ 
eters is given in Tabled 

In order to select “genuine” HAe stars, we applied the 
following selection criteria: 


1. The position in the HR diagram should be in agree¬ 
ment with that of a pre-main-sequence (PMS) star 
with a mass between 1.5 and 3.5 Mq (see Fig.^. 

2. The color criterion J-H>0.25 mag, assuring that the 
onset of the near-infrared (NIR) excess is well defined, 
in agreement with the inner boundary of the gas rich 
disk being set by the silicate evaporation temperature 
of ~1500K (see Table|21for an overview of the photo¬ 
metric data used). 
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Table 2. Literature infrared photometry used in this work. Columns3-7 list the magnitudes in the J (1.25/rm), 
H (1.65/im), K (2.2/rm), L (3.6/im), a nd M (4.8um) pho tomet ric bands, w ith re ferences in colum n 8. T he listed 
reference cod es are : BO: Bouchet et al. |l99ll CA: Carter 1199(1: CO : Cohen ll 97,31: CU: Cutri et al. 1299,31: DW: de 
Wint er et al. l299lt DW2: d e Win ter et al. Il996l: El: Eiroa et al. 299lt FO: Fo uque et al. Il992t GL: Glass & Penston 
Il974t HI: H illenbrand et al. ll992t LA : Lawrence et aljl99^ MA: Malfait et al Jl^Sat ME: Mendoza ll967t ST: Strom 

Columns 9-12 contain the infrared 


et al. Il99(i SY: Sylverster et al. ]l996t VR: Vrba et al. lMld: WA: Waters et al. 


fluxes in Jy from the IRAS Point Source Catalogue I ,Joint IRAS Science Working- GrouDlll988ll . 


(1) 

# 

(2) 

star 

(3) 

J 

(4) 

H 

(5) 

K 

(6) 

L 

(7) 

M 

(8) 

references 

(9) 

12 pm 

(10) 
25 pm 

(11) 
60 pm 

(12) 
100 pm 

1 

AB Aur 

6.10 

5.10 

4.40 

3.30 

2.90 

HI 

27.16 

48.10 

105.60 

114.10 

2 

UX Ori 

8.03 

7.43 

6.71 

5.61 

5.32 

DW,26.2.1986 

2.68 

3.69 

2.85 

~3.76 

3 

HD 36112 

7.44 

6.70 

5.90 

4.75 

4.47 

MA 

5.59 

12.59 

27.98 

18.95 

4 

HK Ori 

9.52 

8.38 

7.29 

5.87 

5.10 

HI 

3.80 

4.08 

<1.64 

<70.37 

5 

HD 245185 

9.34 

8.87 

8.26 

7.46 

6.23 

HI 

4.00 

5.96 

4.97 

3.61 

6 

V380 Ori 

8.25 

7.10 

6.03 

4.44 

3.61 

J-L ST; M ME 

8.61 

8.85 

<75.90 

<38.59 

7 

HD 37357 

8.31 

7.88 

7.24 

6.33 


MA 

2.01 

2.79 

~2.90 

~20.29 

8 

HD 37806 

7.38 

6.68 

5.77 

4.23 

3.79 

MA 

11.02 

9.40 

<5.18 

<33.98 

9 

HD 95881 

7.50 

6.79 

5.87 

4.28 

3.71 

MA 

9.14 

6.87 

1.45 

<1.09 

10 

HD 97048 

7.30 

6.75 

6.04 

4.61 

4.56 

HI 

14.49 

40.34 

69.91 

<250.10 

11 

HD 100453 

6.97 

6.32 

5.52 

4.20 

3.79 

MA 

7.23 

33.59 

39.36 

23.86 

12 

HD 100546 

6.43 

5.88 

5.20 

4.15 

3.80 

MA 

65.78 

242.60 

165.20 

98.56 

13 

HD 101412 

8.70 

8.24 

7.25 

5.81 

5.08 

DW 

3.22 

3.09 

~1.69 

~10.52 

14 

HD 104237 

5.75 

5.14 

4.42 

3.05 

2.58 

MA 

23.65 

23.05 

14.72 

9.58 

15 

HD 135344 

7.44 

6.72 

5.96 

4.76 

4.69 

MA 

1.59 

6.71 

25.61 

25.69 

16 

HD 139614 

7.75 

7.34 

6.76 

5.68 

5.49 

MA 

4.11 

18.14 

19.30 

13.94 

17 

HD 142666 

7.34 

6.72 

6.04 

4.97 

4.69 

MA 

8.57 

11.21 

7.23 

5.46 

18 

HD 142527 

6.65 

5.94 

5.20 

3.89 

3.50 

MA 

10.38 

21.23 

105.10 

84.70 

19 

HD 144432 

7.21 

6.69 

6.14 

5.14 

4.90 

MA 

7.53 

9.36 

5.76 

3.29 

20 

HD 144668 

5.83 

5.18 

4.38 

3.08 

2.54 

HI 

18.05 

14.51 

~14.36 

<63.25 

21 

HD 150193 

7.05 

6.37 

5.64 

4.37 

3.93 

MA 

17.61 

18.10 

8.13 

<16.25 

22 

HD 163296 

6.24 

5.52 

4.70 

3.52 

3.14 

MA 

18.20 

20.99 

28.24 

<40.62 

23 

HD 169142 

7.43 

7.01 

6.53 

5.64 

5.57 

MA 

2.95 

18.43 

29.57 

23.42 

24 

HD 179218 

6.99 

6.64 

5.91 

4.68 

4.18 

J,H El; K-M LA 

23.44 

43.63 

29.92 

17.35 


2.1. Selection effects 

In order to identify if other selection effects may have 
entered our sample, we determined the mass and age 
(calculated from the birthline, for an accretion rate of 
19“®M0yr“^) of the stars by placing them in the HR 
diagram (Fig.P) and comparing their positions to PM S 
evolutionary tracks published by IPalla fc Stahlej lll99,3ll . 
Throughout this work “(PMS) age” r is defined as the 
time that past since the star was on the birthline, which 
is when the star becomes optically visible for the first 
time. The luminosity is calculated from the observed pho¬ 
tometry and the measured distance to each star. The un¬ 
certainties in the derived stellar masses and ages listed 
in Tabled] and shown in Fig. [21 mainly reflect the un¬ 
certa inty in the distances (see also lyan den Ancker et af] 
Il99j1l . A comparison of the st ellar parameters deriyed us¬ 
ing the lPalla fc Stahl^ lll99,3^ evolutionar y tracks to those 
obtained using more r ecent calculations llPalla fc Stahleil 
I 1999 I: ISiess et ^ 129901) ■ shows that the particular choice 
of PMS tracks introduces an additional uncertainty in the 
stellar ages that is of the same order as those implied 
by observational uncertainties in the position of the star 


in the HR diagram. This additional uncertainty is not 
taken into account. In Fig.|21 we show the derived stel¬ 
lar masses versus the PMS ages of the stars. It is evident 
from this figure that our sample lacks ’’old” (r > 19®'^yr) 
disks around 3-4Mq stars, and ’’young” (r < 19®yr) disks 
around stars less massive than about 2.5 Mq. The lack of 
“old” disks around 3-4 Mq stars is likely caused by the 
shorter timescale on which the disks around more massive 
stars are dispersed. It is therefore not a selection effect. 
The lack of “young” disks around the lower mass stars is 
most likely caused by the fact that these stars do not clear 
their environment as rapidly as the more massive (more 
luminous) stars, and consequently become optically bright 
later in their evolution. Also, at a given age, lower mass 
PMS stars are less luminous than more massive ones, and 
so more easily escape optical detection. This selection ef¬ 
fect must be taken into account in any discussion about 
the evolution of the dust in proto-planetary disks based 
on optically selected samples. 
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Fig. 1. The positions of our stars in the HR diagram. The 
solid curves indicate the pre-main sequence evolutionary 
tracks of stars of different masses; the dashed curve repre¬ 
sents the birthline fo r an accretion rate of 10 “® M 0 yr“^ 
(both are taken from IPalla fc Stahleilll99^ . 


4.0 


3.5 


3.0 


G 


° 2.5 


2.0 


1.5 


1 r- 


I I 


1.0 I_,_u 


5.5 



6.0 6.5 7.0 

log(PMS age) [log yr] 


7.5 


Fig. 2. The derived stellar mass estimates vs. the pre- 
main-sequence age estimates derived from a comparison 
of their position in the HR diagram to PMS evolutionary 
tracks of Palla & Stabler (1993). Our sample shows a clear 
lack of disks around relatively old stars of 3-4 Mq, and 
of disks around relatively young stars less massive than 
about 2.5 Mq. 

2.2. Classification of the sources 

MEOl empirically decomposed the infrared spectra of 
Herbig Ae/Be stars into three components: a power law 
component, a cold black-body component, and solid state 
emission bands (mainly at 10 and 20 fim). They found that 
some sources exhibit both the power law and cold black- 
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_l 

Z 

_l 

1 


Fig. 3. Classification of the sources based on global SED 
properties. We plot the ratio of the near-infrared and in¬ 
frared luminosity (see section vs. the IRAS mi 2 -m 6 o 
color (defined as mi 2 -m 6 o=—2.5logiFi2/.7^60, where T 12 
and J-Qo are the fluxes at 12 and 60 /rm as listed in the 
IRAS point source catalogue). Sources classified as group I 
according to MEOl are in the lower right part of the dia¬ 
gram, group H sources are in the upper left part. Group la 
sources are indicated with crosses, group Ib sources by tri¬ 
angles, and group H sources with diamonds. The dashed 
curve indicates our division line between the two groups 
('section l2.2ll . 


' ' ' ' I. I .I.I. I 

_ group II 


o20 



I. I.I .I 


- 2-10123 

^ 12“^60 


body component, and classified these sources as “groupI”. 
Sources that lack the cold black-body component were 
classified as “group H”. A further division into subgroups 
“a” and “b” serves to indicate the presence or absence, 
respectively, of silicate emission bands at 10 and 20 fj.ia. 
It was proposed by MEOl that group I sources have a 
large (several hundred AU) flared outer disk, whereas the 
group H sources have a smaller, non-flaring outer disk. 

We classify the sources for which we have newly 
measured N-band (8-13.5p,m) spectra following MEOl. 
Whereas MEOl had ISO spectra of their sources at 
their disposal, our classification is based solely on broad¬ 
band photometry. We find that the group I and group H 
sources are well separated in an IRAS mi 2 -m 6 o color ver¬ 
sus Lnir/Lir diagram (Fig. EJ, where Lnir is the inte¬ 
grated luminosity as derived from the J,H,K,L and M 
band photometry, and Lir is the corresponding quan¬ 
tity derived from the IRAS 12, 25 and 60/rm points. 
For groupl sources, Lnir/Lir < (mi2-m6o)-l-1.5, groupll 
sources have Lnir/Lir > (mi2-m6o)-l-1.5. The dashed line 
in Fig. 13 indicates the boundary between the two groups. 
To avoid any pre-assumptions concerning the spectral 
shapes of our targets, we did not apply color correction 
to the IRAS data for the classific ation of the sources. 

By use of the Lnir/Lir ratio llvan Boekel et al ] |2003ll 
we compare the near-IR SED to the mid-IR SED. The 
near-IR SED is found to be similar for all HAe stars 
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llNatta et alJl200ili . while the major differences occur in 
the mid-IR SED. The TnirZ-Z^ir ratio is smaller for group I 
than for group II sources. The mid-IR SED of group I 
sources is “double-peaked” compared to the SED of a 
group II source. Group I sources are redder than their 
group II counterparts. The IRAS mi 2 -m 6 o color index 
quantifies this difference in SED shape. 


3. Observations and data reduction 

Infrared spectra in the 10 /rm atmospheric window were 
taken in December 2001 and March 2003 with the 
TIMMI2 instrument mounted at the 3.60 m telescope at 
ESO La Silla observatory. Conditions were clear during 
all nights. The low resolution (RwlGO) N band grism was 
used in combination with a 1.2 arcsecond slit. The pixel 
scale in the spectroscopic mode of TIMMI2 is 0.45 arcsec- 
onds. A log of our observations is given in TablelHl 

3.1. Atmospheric correction 

Ground-based observations at 10 pm are particularly chal¬ 
lenging because of the high atmospheric and instrumental 
background, and the varying transmission of the earth at¬ 
mosphere. As proper spectral calibration is essential for 
the compositional analysis presented in section IKdl we re¬ 
port here in some detail on our data reduction method. 

To deal with the high background signal, we employed 
standard chopping and nodding, using a -1-10 arcsecond 
chop throw North-South, and a —10 arcsecond nod throw 
North-South. For the spectral calibration of our measure¬ 
ments we regularly observed standard stars. These obser¬ 
vations are used to determine both the atmospheric ex¬ 
tinction per unit airmass (A^) and the instrumental re¬ 
sponse (i?y) at all wavelengths. For each science observa¬ 
tion, we determine Ai, and R^, from two calibration mea¬ 
surements: 


Su.l — li^.l ^ ^1^,1 — Ai^TTZA,! (1) 

Siy.2 — C ’ Ri/, ^1^,2 — Aiy?7ZA,25 

where Sj^j are the measured calibrator spectra, Rj are 
the intrinsic (model) calibrator spectra, are the opti¬ 
cal depths of the earth atmosphere during the calibration 
measurements, and rriA.i denote the airmass at which the 
calibrators were observed. Whenever possible, we chose 
rriA,! or toa ,2 to be very close to the airmass of the sci¬ 
ence observation. Both calibrator measurements are cho¬ 
sen as close in time as possible to the science observation. 
For the in trinsic calibrator sp ectra we use “spectral tem¬ 
plates” bv iGohen et al.l lll999ll ^. Tabled lists the observed 
calibrators and the applied templates, which are chosen to 

^ available for download at the TIMMI2 website 
(http://www.Is.eso.org/lasilla/sciops/timmi/docs/tables/ ) 


match the calibrator spectral type as closely as possible. 
Solving Eq. I for A^, and R^,, we find: 


A. 


ln(R, 2 /R,i) + ln(s;.,i/si., 2 ) 

toa ,2 - rriA.i 


( 2 ) 


Rv 






— AifiTiA ; 


(3) 


The intrinsic spectrum of a science target observed at air¬ 
mass rriA is then calculated from its measured spectrum 
Sy as: 


T _ ^ ^Ai/TTlrmA 

lu — e 
ixi> 


(4) 


For a number of sources we have ISO spectra available. 
We generally find very good agreement in spectral shape 
between the ISO data and our new ground based spectra. 
A comparison with the ISO spectra of the brighter stars 
(HD 100546, HD 163296, HD 150193) shows differences at 
the level of at most a few percent in the shape of the spec¬ 
tra, when the TIMMI2 spectra are scaled such as to most 
closely match the flux levels in the ISO spectra. For each 
HAe star spectrum, the time and airmass of observation, 
as well as the time and airmass of the two calibration ob¬ 
servations are given in Table|31 


3.2. Flux calibration 

Our stars have been selected to be isolated Herbig stars. 
Therefore, the IRAS photometry will in most cases not 
be contaminated by emission from nearby sources or 
a surrounding remnant cloud. All emission seen in the 
IRAS data is expected to originate in the disk. We flux- 
calibrated the spectra using the IRAS 12 ^m data, by ap¬ 
plying a scaling factor q such that 

^OO 

q / = iFi 2 X 1.35 X 10~^^ [erg s“^cm“^], (5) 

where Ti, is the normalized instrumental response func¬ 
tion of the IRAS 12 micron band, and iF \2 is the 12 mi¬ 
cron flux as listed in the IRAS point source catalogue (a 
source that has a 12 pm flux of J -\2 = l Jy in the IRAS point 
source catalogue, yields an inband flux of 1.35 x 10“^° 
erg s“^cm“^). The IRAS 12 micron band runs from 8 to 
15 ^m and is therefore somewhat broader than our spec¬ 
tral coverage. To allow for the calibration we estimate the 
spectrum between 13.5 and 15 pm, by linearly extrapolat¬ 
ing R between the continuum points measured at 8 and 
13 ^m. Our final calibrated spectrum is then E/ = qR. We 
estimate the absolute photometric accuracy to be 15%. 


4. Description of the observations 

4.1. Description of the spectra 

The spectra of our sample of HAe stars are shown in 
Fig.EI All stars show spectral structure on top of a con¬ 
tinuum whose slope varies strongly from star to star. 
The spectral features are due to various kinds of silicates 
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Table 3. Log of the TIMMI2 observations. We list the observing date (defined as the day on which each observing 
night began), time (U.T.), airmass of the observation, and integration time in seconds (columns 3-6). The calibrators 
used for the atmospheric correction are also given, with the time and airmass of the measurements (columns 7-12). 


(1) 

# 

(2) 

star 

(3) 

date 

(4) 

time 

(5) 

rriA 

(6) 

Tint 

(7) 

calibrator 1 

(8) 

time 

(9) 

mA,i 

(10) 

calibrator 2 

(11) 

time 

(12) 

mA,2 

2 

UX Ori 

27-12-2001 

01:09 

1.31 

1380 

HD 32887 

01:59 

1.06 

HD 32887 

00:39 

1.25 

3 

HD 36112 

18-03-2003 

23:19 

1.83 

960 

HD 48915 

00:27 

1.04 

HD 29139 

23:58 

1.88 

4 

HK Ori 

19-03-2003 

23:21 

1.39 

960 

HD 48915 

00:56 

1.08 

HD 23249 

00:04 

1.68 

5 

HD 245185 

19-03-2003 

00:27 

1.56 

720 

HD 48915 

00:56 

1.08 

HD 23249 

00:04 

1.68 

6 

V380 Ori 

17-03-2003 

00:59 

1.37 

960 

HD 48915 

00:21 

1.04 

HD 58972 

01:36 

1.34 

7 

HD 37357 

19-03-2003 

01:08 

1.43 

720 

HD 48915 

00:56 

1.08 

HD 55865 

02:33 

1.46 

8 

HD 37806 

27-12-2001 

04:04 

1.12 

690 

HD 32887 

01:59 

1.06 

HD 32887 

00:39 

1.25 

9 

HD 95881 

17-03-2003 

02:55 

1.37 

720 

HD 48915 

00:21 

1.04 

HD 107446 

02:41 

1.32 

10 

HD 97048 

19-03-2003 

03:02 

1.52 

720 

HD 98292 

04:02 

1.28 

HD 55865 

02:33 

1.46 

11 

HD 100453 

18-03-2003 

05:24 

1.13 

960 

HD 93813 

05:00 

1.08 

HD 146003 

06:38 

1.25 

12 

HD 100546 

19-03-2003 

05:09 

1.34 

960 

HD 109379 

05:36 

1.01 

HD 98292 

06:18 

1.37 

13 

HD 101412 

17-03-2003 

03:25 

1.19 

960 

HD 48915 

00:21 

1.04 

HD 89388 

04:04 

1.19 

14 

HD 104237 

17-03-2003 

04:57 

1.52 

480 

HD 123139 

08:04 

1.03 

HD 89484 

04:37 

1.65 

15 

HD 135344 

17-03-2003 

05:56 

1.13 

960 

HD 123139 

08:04 

1.03 

HD 107446 

06:35 

1.20 

16 

HD 139614 

17-03-2003 

07:25 

1.05 

960 

HD 123139 

08:04 

1.03 

HD 107446 

06:35 

1.20 

17 

HD 142666 

18-03-2003 

06:07 

1.23 

720 

HD 93813 

05:00 

1.08 

HD 146003 

06:38 

1.25 

18 

HD 142527 

18-03-2003 

06:56 

1.11 

720 

HD 123139 

08:07 

1.04 

HD 146003 

06:38 

1.25 

19 

HD 144432 

18-03-2003 

07:25 

1.06 

960 

HD 123139 

08:07 

1.04 

HD 146003 

06:38 

1.25 

20 

HD 144668 

17-03-2003 

09:15 

1.02 

720 

HD 152334 

09:45 

1.03 

HD 152786 

09:01 

1.14 

21 

HD 150193 

18-03-2003 

08:23 

1.04 

720 

HD 123139 

08:07 

1.04 

HD 146003 

06:38 

1.25 

22 

HD 163296 

18-03-2003 

09:40 

1.04 

480 

HD 123139 

08:07 

1.04 

HD 152786 

09:22 

1.12 

23 

HD 169142 

17-03-2003 

10:07 

1.03 

960 

HD 152334 

09:45 

1.03 

HD 152786 

09:01 

1.14 

24 

HD 179218 

19-03-2003 

09:39 

1.78 

480 

HD 152334 

07:59 

1.10 

HD 187642 

10:00 

1.68 


Table 4. Spectral templates used for the calibrators (see 
section lO) . 


calibrator 

spectral 

type 

template file 

template 
sp. type 

HD 23249 

KOIV 

hdl23139.tem 

KOIIIb 

HD 29139 

K5HI 

alp_Tau.dat 

K5HI 

HD 48915 

AIV 

law 

law 

HD 55865 

KOHI 

hdl23139.tem 

KOIIIb 

HD 58972 

K3HI 

hd6805.tem 

K2HI 

HD 89388 

K3IIa 

hd6805.tem 

K2HI 

HD 89484 

Kllllb 

hdl69916.tem 

KlHIb 

HD 92397 

K4.5HI 

hd32887.tem 

K4HI 

HD 93813 

KO/KIIH 

hdl23139.tem 

KOIIIb 

HD 98292 

M2HI 

alp_Tau.dat 

K5HI 

HD 107446 

K3.5HI 

hd32887.tem 

K4HI 

HD 109379 

G5H 

hd37160.tem 

G8HI-IV 

HD 123139 

KOIIIb 

hdl23139.tem 

KOIIIb 

HD 146003 

M2HI 

alp_Tau.dat 

K5HI 

HD 152334 

K4HI 

hd32887.tem 

K4HI 

HD 152786 

K3HI2 

hd6805.tem 

K2HI 

HD 177716 

KlHIb 

hdl69916.tem 

KlHIb 

HD 187642 

A7V 

HD187642 spec.dat 

A7IV-V 


(see section 15.11 below). Also, emission from Polycyclic 
Aromatic Hydrocarbons (PAHs) can be seen, at 7.9, 8.6, 
11.3 and 12.7/xm. The relative importance of the silicate 


and PAH contributions varies strongly. There are sources 
that show both silicate and PAH emission (e.g. HD 100546, 
HD 179218), sources that show only silicate emission (e.g. 
HD 144432), and sources that display only PAH emis¬ 
sion (these are the group Ib sources HD 97048, HD 100453, 
HD 135344, and HD 169142). There are no sources in our 
sample that have a completely featureless 10 micron spec¬ 
trum. 

MEOl and lAcke fc van den Anckeil ll2004^ found that 
the SED correlates with the presence and/or strength of 
the PAH bands: group 1 sources tend to show (prominent) 
PAH emission, while group H sources do not. This trend 
is confirmed in our sample, but we note that there is con¬ 
siderable scatter. For instance, HD 95881 has little far-IR 
excess and is thus classified as group 11, but, nonetheless, 
shows clear PAH emission bands. 

The silicate band shows very large variations in shape 
and strength. The bulk of the emission is in most cases 
due to amorphous silicates, but almost all stars show some 
spectral structure near 11.3/xm, which can be attributed 
to forsterite. Note, however, that this feature blends with 
the 11.3/xm PAH band. There are also prominent narrow 
emission bands near 9.2 and 10.6/xm. These are due to 
crystalline enstatite. The spectrum of HD 100546 is domi¬ 
nated by crystalline forsterite, while that of HD 179218 is 
dominated by crystalline enstatite. This latter star shows 
one of the richest 10 ^m spectra observed to date (see 
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FigEJ. The resonances of crystalline enstatite are clearly 
visible in the spectrum of this source. The ISO spectrum of 
HD 179218 at longer wavelengths also poi nts to a relatively 
high abundance of crystalline enstatite I Bouwman et al.| 
l2nmli . It is obvious that the nature of the crystalline dust 
in our sample shows very large variations, both in terms 
of the fraction of the dust that is crystalline, and in com¬ 
position. 

There are four stars (HD 97048, HD 100453, HD 135344 
and HD 169142), all classified as group I, that show no de¬ 
tectable silicate emission. Instead, their 10 /rm spectra are 
dominated by PAH emission. The lack of silicate emission 
is most simply explained by assuming that there are no 
small (< 3-5 »m) silicate gr ains in the inner 10-20 AU of 
the disk. iMeeiis et al.l ( 2002^1 derive limits on the presence 
of small silicate grains in HD 100453, and argue that all 
grains smaller than 4/rm must have been removed. The 
most likely cause for the removal of small silicate grains is 
grain growth, but apparently this has not affected the pop- 
ulation of small carbona ceous grains to the same extent 
llvan Boekel et al.l2004lJl . Possibly, small grains survive in 
the outer disk regions. At large distance from the star, the 
silicate grains may be too cold to contribute to the 10 /rm 
spectrum, while the PAHs can still produce significant 
emission. This can only occur if the PA Hs have not been 
incorpo rated into larger grains. Indeed, Ivan Boekel et al.l 
ll2nn4ah find evidence for a distance dependence of the 
typical silicate grain size in the surface layers of HAe star 
disks: in the innermost regions growth has proceeded fur¬ 
ther than in the outer disk regions. 


4.2. The shape and strength of the silicate feature 


In Fig.l^we show the ratio of the normalized spectrum at 
11.3 and 9.8 pm, against the silicate peak to continuum 
ratio (“feature strength”)^. Group I sources are indicated 
with triangles, group H sources with diamonds (the “out¬ 
lier” ff 12 is HD 100546). Sources with a low 11.3/9.8 ratio 
have a triangular shaped emission feature, clearly peaked 
just shortward of 10 ^m (e.g. UX Ori, in Fig. 0 . Sources 
with a high 11.3/9.8 generally have a broad, flat-topped 
emission band, often showing substructure (e.g. HD 37806, 
HD 142527). There is a clear correlation between the shape 
and strength of the silicate feature; stars with a strong fea¬ 
ture (i.e. a high peak/continuum ratio) have a low 11.3/9.8 
ratio, whereas stars with weaker silicate features have 
higher 11.3/9.8 ratios. This cor relation was first demon ¬ 


strated in Herbig Ae stars by TOi^^oete^_e^^ ] (H^. 


iMeeiis et al ] (I2OO3II and IPrzvgodda et alJ ll 200 ,‘ji ' have sub- 


^ To estimate the continuum we simply interpolate linearly 
between 8 and 13 micron. The peak/continuum ratio is the 
maximum value of the normalized spectrum Fnorm = 1 -|- 
Fv,cs/ < Fi,^c >, where is the continuum subtracted spec¬ 
trum (G/ — F,/,c) and < G/,c > is the mean of the continuum. 
This definition of Fnorm preserves the shape of the emission 
band even if the continuum is not constant. For a constant 
continuum level it is identical to Fv/F^.o. 



Fig. 5. The N-band spectrum of HD 179218 (upper 
panel), and the measured mass absorption coeffic ients of 
ortho enstatite taken from IChiha.ra, et alJ ( 200211 (lower 
panel). The wavelengths of the most prominent emission 
bands are indicated by the dotted lines. In this object, 
enstatite grains are an important constituent of the grain 
population that causes the 10 pm feature. 


sequently shown that the same trend is observed in the 
silicate feature of T-Tauri stars. 

The shape of the emission bands with a low 11.3/9.8 
ratio is similar to that of the ISM silicate absorption 
feature. Such 10 micron features are indicative of small, 
amorphous silicate grains, i.e. relatively “primitive” dust. 
The emission bands with high 11.3/9.8 ratios can be ex¬ 
plained with on average larger grains, and a higher de¬ 
gree of crystallinity, i.e. relatively “processed” dust. Thus, 
the silicate feature 11.3/9.8 ratio is a measure of the 
amount of processing that the material has undergone 
I Bouwman et aLlboOlll . 

5. Analysis 

5.1. Compositional fits 

To derive the composition of the silicate dust causing 
the 10 micron feature, the observed spectra have been 
fitted using the most commonly found dust species in 
circumstellar material that show spectral structure in 
the 10 micron region (see Fig. [7|). These are amorphous 
and crystal line olivine and pyrox ene, and amorphous 
silica (e.g. iBouwman et al ] 1200111. Amorphous olivine 
(Mg 2 a;Fe 2 (i_a;)Si 04 , where 0 < x < 1 denotes the mag¬ 
nesium content of the material) is the most commonly 
found silicate in astrophysical environments. It dominates 
the 10 micron extinct ion caused by dust grains in the ISM 
I Kemper et al ] |2004 . The 10 micron emission spectrum 
of small amorphous olivine grains is characterized by a 
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Fig. 4. N-band spectr a of the source s in ou r sample. The ISM silicate extinction efUcie ncv, plotted in the upper left 
panel, was taken from iKemner et ^ ll2flfl4l . The AB Aur spectrum was taken by ISO I van den Ancker et alJEflfinjl . 
Also plotted are the best fits to the spectra (grey curves, see section lb^ . 


rather broad feature which peaks at 9.8 /rm. Small amor¬ 
phous pyroxene grains (MgajFei-ajSiOa) show an emission 
feature very similar to that of amorphous olivine grains, 
though shifted toward shorter wavelengths. The emis¬ 
sion spectra from small crystalline olivine and pyroxene 
grains show strong, narrow resonances that are observed 


1998 

Bouwman et al. 

1222 : 

Bouwma.n et al. 



and comets |Crovisie7'end 


l200^ . From the positions of the 


resonances in the emission spectra of these objects it is 


cnesium rich 

compt 

l.lae-er et all 



In our fitting pro¬ 
cedure we therefore use crystalline olivine and pyroxene 
with a: = I, i.e. forsterite and enstatite respectively. For 
the amorphous olivine and pyroxene we use x = 0.5. 

_Studies of interplanetary dust particles I R ietmeii^ 
have shown that some of these particles contain 


large inclusions of silica (Si02). Also, from laboratory ex¬ 
periments it is suggested that when amorphous silicates 
are annealed to form fo rsterite, silica will be formed (see 
e.g. iFabian et alJl200nll . The emission spectrum of silica 
in the 10 micron region has a distinct spectral signature, 
with a strong feature peaking at 8.9/rm. Therefore, silica 
is included as one of the possible dust components. 

Many of our sources show emission bands at 7.9, 8.6, 
11.3 and 12.7/xm that are attributed to PAHs. In order 
to include the PAH emission in our compositional fits, we 
constructed a simple PAH template. This was done by 
taking our two highest quality spectra of sources without 
a silicate feature (HD 97048 and HD 169142), subtracting 
the continuum emission, and averaging over the two spec¬ 
tra. The resulting PAH emission spectrum was added as a 
fit component, and is shown in the bottom panel of Fig.^l 
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Fig. 4. (Continued) 


Dust grains in circumstellar environments are most 
likely aggregates of many different dust species. However, 
in order to perform an analysis of the abundances of the 
various components in a large sample of sources, it is not 
(yet) feasible to do computations for such complex aggre¬ 
gates. If the aggregates are very fluffy, the constituents 
of which they are composed will i nteract with the inci- 
dent radiation as separate entities. iMolster et alJ ll200,‘lll 
show that indeed a measurement of the infrared spec¬ 
trum of a relatively large IDP still displays the spectral 
structure one would expect from much smaller particles. 
Therefore, we assume that the emission properties of an 
aggregated structure can be represented by the sum of the 
emission prop erties of the constituents (for a similar ap¬ 
proach see e.g.lfeucato et ahlllQQllI: iBouwman et al.ll200lt 
iHonda et akfetin^ 

All dust grains contributing to the 10 ^m region will 
in principle have their own temperature depending on the 
grain size and composition as well as on position in the 
disk. Consequently, the emission function of the dust will 
be a (weighted) sum over blackbody curves of different 
temperatures. In order to account for these effects one 
should use a disk model that self-consistently describes the 
disk geometry and all relevant radiation processes. This is 
beyond the scope of the analysis p resented here. For state 
of the art disk models we refer to IPiiHemond DominilJ 


il2nn4 . We add that even these sophisticated models still 
do not account for many relevant processes, including the 
settling and (turbulent) mixing of grains leading to spa¬ 
tial gradients in the d ust properties. Such spatia l gradients 
have been reported bv ivan Boekel et ^ ( 2004a^ to exist in 
Herbig Ae disks. One should also realize that even if such 
a detailed modeling approach was attempted, the modest 
wavelength interval provided by the TIMMI2 data would 
reveal only very limited information on the temperature 
distribution. The continuum in the 10 micron region is 
a very smooth, almost linear function of wavelength. We 
therefore opted to represent the continuum by a single 
blackbody curve with a temperature, that is charac¬ 
teristic for the dust emission as a whole. Inherent in this 
approach is that we essentially assume that the tempera¬ 
ture of the individual grains is independent of composition 
and size, and that there are no gradients in grain proper¬ 
ties (per unit mass) throughout the disk. 

The disk regions that we study are partially optically 
thick at 10 ^m. However, we assume that the disk surface 
layer from which the observed flux originates is optically 
thin at this wavelength. This must be the case since we 
see emission features. The flux emitted by a distribution 
of dust grains is then given by 




silicate 


cx B,{T,) V Win 
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Peek over Continuum 

Fig. 6. The flux ratio of the normalized spectra at 11.3 
and 9.8 micron (a measure for the amount of pro¬ 
cessing that the material has undergone) versus the 
peak/continuum ratio of the silicate feature (a measure 
for the typical grain size). Group I sources are represented 
by triangles, groupII sources by diamonds. In the upper 
right corner of the hgure we have indicated the typical 
uncertainties in the displayed quantities. 


where B^iTc) denotes the Planck function at the charac¬ 
teristic temperature Tc, Ki is the mass absorption coeffi¬ 
cient of dust component i (see also Fig. 0 ), and Wi is a 
weighting factor which is proportional to the total dust 
mass in component i. The summation is over all dust 
components. The mass absorption coefficient of each dust 
species is determined by the size, shape, structure and 
chemical composition of the dust grains. The total model 
spectrum jg calculated by adding a continuum 

and PAH contribution to the silicate emission. 

5.1.1. Shape of the dust grains 

The shape and structure of the dust grains are very im¬ 
portant parameters determining the feature shape of the 
emission spectrum. Usually it is assumed that the grains 
are homogeneous and spherical so that Mie theory can 
be applied to calculate the Ki. Another widely used as¬ 
sumption is that the grains are much smaller than the 
wavelength of radiation (“Rayleigh limit”) in which case 
it is mathematically straightforward to adopt a contin¬ 
uous distribution of ellipsoi ds (CDE iBohren fc Hiiffm'ai] 
Il983l iBouwman et alJ^2003^ . Since micron sized silicate 
grains are not in the Rayleigh limit at a wavelength of 
10 /im, we cannot use CDE calculations to study grain 
growth. Furthermore, a comparison between calculations 
of the mass absorption coefficients of small crystalline sil¬ 
icates with measurements shows that we cannot get good 


agreeme nt using homogeneo us spherical particles (see for 
example iFabian et al.ll200ll) . 

Adopting different grain shapes, iMin et akl ll2003ll 
showed that the absorption properties can be divided in 
essentially two categories. One category contains the per¬ 
fect homogeneous spheres; the other all other investigated 
shapes, including hollow spheres. Shape effects within the 
second category do exist, but they are small compared to 
the differences with homogeneous spheres. One could say 
that the difference between perfect homogeneous particles 
and those having other shapes is essent ially a result of a 
breaking of perfect symmetry (see also iMin et 
This implies that we have only very limited information 
on the true (likely irregular) shape of astrophysical dust 
grains from spectroscopic analysis. However, a practical 
implication of this result is that one may represent the ab¬ 
sorption properties of irregular grains with sufficient accu¬ 
racy by adopting the average properties of a distribution 
of shapes other than that of homogeneous spheres. For 
this purpose, a practical choice is a distribution of hollow 
spheres (DHS), simply averaging over the volume fraction 
occupied by the central inclusion, which ranges from 0 to 
1. In this shape distribution, the material volume of the 
particle is kept constant, thus particles with a high value 
of / will have a large outer radius. This shape distribu¬ 
tion has the advantage that it can be applied f or all grain 
sizes u sing a simple extension of Mie theory. iMin et ffil 
1 I 2 OO 3 I) showed that this indeed gives excellent results for 
small forsterite grains. Therefore, in this work, the mass 
absorption coefficients of all crystalline grains (forsterite, 
enstatite) and silica are calculated with a distribution of 
hollow spheres. For the amorphous olivine and pyroxene 
particles we use Mie theory since for these species the ef¬ 
fects of shape on absorption properties are minor. 

The Ki are calculated using laboratory measurements 
of the refractive index as a function of wavelength. 
References for the measurements used for the various dust 
species are listed in TableO 

5.1.2. Size of the dust grains 

The dust grains in circumstellar disks most likely have a 
rather broad size distribution. In the 10 micron region 
the observational data are sensitive to the dust grains 
with a volume equivalent radius up to a few micron in 
size. Larger grains mainly contribute to the continuum. 
In order to minimize the number of free parameters in 
the htting procedure, we want to sample the size dis¬ 
tribution carefully. We hnd that the variety of spectral 
shapes can be best covered using only two distinct par¬ 
ticle sizes, a ‘small’ particle size with a material volume 
equivalent sphere radius ry = 0.1/rm and a ‘large’ par- 
ticle size with ry = 1 .5 qm (for a simila r approach see 
iBouwman et al.1 [20011 iHonda et al.ll2003i . We have ex¬ 
tensively checked that the results of the analysis using 
three, four or five particle sizes with volume equivalent 
radii ranging up to 3.5/im do not change significantly. 
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Dust component 

Chemical formula 

Shape 

Reference 

Amorphous olivine 
Amorphous pyroxene 
Crystalline forsterite 
Crystalline enstatite 
Amorphous silica 

Mg2i:Fe2-22:Si04 

Mg,cFei_a:Si03 

Mg2Si04 

MgSiOa 

Si02 

Homogeneous Spheres 
Homogeneous Spheres 
Irregular (DHS) 
Irregular (DHS) 
Irregular (DHS) 

Dorschner et al^9^ 
Dorschner et al^9^ 
Servoin & Piriou^^^ 
Jager et al.^ffi^ 

Spitzer & Kleinman^9^ 


Table 5. Characteristics of the various dust components used in the fitting procedure (see section |OJ. The chemical 
formulae and assumed grain shapes have been indicated in the second and third column. In the fourth column, we 
give references for the refractive index data used. 


The size of the large grain component, when using only 
two particle sizes, has to be chosen with care. When it 
is too large, the difference in absolute value of the mass 
absorption coefficients of the small and the large grain 
component will be too big. If, for example, one adopts 
for the large grain size ry = 2 /rm, one would overesti¬ 
mate the abundance of large grains. This is especially so 
for the large crystalline grains and results in an overes¬ 
timate of the crystallinity. The differences between the 
absolute values of the emissivities of a 0.1 and a 2/rm 
amorphous grain are not so big. Actually, adopting 2 /rm 
sized grains is even slightly better for the l arge amorphous 
component (see also lBouwman et al .12001 h . However, con¬ 
sidering both the amorphous and crystalline component, 
1.5 ^m is the best choice for the large grains. 

5.1.3. Fitting procedure and error analysis 

In order to keep the number of free parameters in the 
model small, only two grain sizes are used for every sili¬ 
cate dust type (as discussed above). We thus have contri¬ 
butions of five silicate species, of PAHs and of a blackbody 
continuum of which the absolute level and shape (charac¬ 
teristic temperature) can be varied. The emissivities of the 
silicates are multiplied by a blackbody spectrum with the 
same characteristic temperature as the continuum. This 
results in 13 free parameters. The silicate and PAH tem¬ 
plates are shown in Fig.H 

To fit the spectra we minimize the reduced of the 
entire 10 micron region given by 

.. 'T'model f \ \ ^observed / \ ^ 

2 _ ^ '’u _ /r,\ 

^ ~ Nx-M^ a, ■ ^ ^ 

Here N\ is the number of wavelength points Ai, M is the 
number of fit parameters (in this case M = 13) and ai is 
the absolute error on the observed flux at wavelength Ay 
For a given characteristic temperature we can calculate 
the optimal values for the weights Wi of the individual dust 
components, using a linear least square fitting procedure. 

The measurement errors (ct^) used in the fitting pro¬ 
cedure represent the statistical noise in the spectra. The 
calibrator spectra all have a very high signal to noise ratio 
(SNR) and statistical noise of the calibration observations 
is negligible. The SNR in our Herbig star spectra range 
from ^18 in the faintest source (HD 135344) to approxi¬ 
mately 60 in the bright sources (e.g. HD 100546). There 


are also systematic uncertainties, arising from an imper¬ 
fect calibration, and the uncertainty in the used spectral 
templates for the calibrators. Some degree of systematic 
error is inevitable, since the science target and calibrator 
cannot be measured at the same time and in the same 
direction. Since we cannot assess the systematic uncer¬ 
tainties we do not take these into account. We note, how¬ 
ever, that agreement between our ground based spectra 
and high SNR ISO spectra is generally very good. For the 
faint sources (such as UX Ori), the statistical noise domi¬ 
nates the error budget. For bright sources (e.g. HD 144432) 
the systematic uncertainties may be important, implying 
that we underestimate the errors. This will evidently lead 
to higher values in the fit procedure. 

The errors on the fit parameters are calculated using 
a Monte Carlo method. For every spectrum we gener¬ 
ate 1000 synthetic spectra, by randomly adding Gaussian 
noise to the spectrum with a distribution of width ai at 
each wavelength point. This yields 1000 spectra that are 
all consistent with our data. On each of these, we per¬ 
form the exact same compositional fit procedure, yielding 
(slightly) different values of the fit parameters. From the 
resulting distribution of all fit parameters, we calculate 
the mean (which will be our ’best fit’ value) and stan¬ 
dard deviation. Besides its simplicity, this method has the 
advantage that degeneracies between fit parameters auto¬ 
matically show up as large errors in these parameters. 

5.2. Results 

Fig.m shows the observed spectra together with the best 
fit model spectra. For comparison we also fit the inter- 
stellar extinction as observed towards the galactic center 
llKemoer et ahlEoP^ . The resulting values of the fit pa¬ 
rameters are summarized in TableEl The overall quality 
of the model fits is very good and we thus conclude that 
the diversity of the shapes of the observed spectra are well 
covered by the choice of the spectral templates. We notice 
that in some sources the model fits show slightly more 
spectral structure between 9 and 11 fim than the observed 
spectra. This could be caused by the choice of the shape 
distribution of the dust grains, which may be too sim¬ 
ple to represent the spectral details of realistic particles. 
We also notice that in, for example, HD 179218 we can¬ 
not accurately reproduce the detailed shape of the spec¬ 
trum. Especially the feature around 11.3 /im is less sharply 
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A [ytim] 

Fig. 7. The mass absorption coefficients of the various 
templates used in the fitting procedure (upper 5 panels). 
We use grains with volume equivalent radii of 0.1 /rm (solid 
lines) and 1.5/im (dotted lines). In the lower panel we 
show the template used for the PAH emission, which is 
normalized such that the maximum flux in the 8 to 13 
micron region equals unity. For a detailed discussion see 
section o 


peaked in the observed spectrum than in the model fit. 
This effect can also be seen in other sources, albeit in a 
more modest form, and could be caused by a missing dust 
component or by the choice of the shape distribution. 

As an objective measure for the goodness of fit, the re¬ 
duced of every fit is listed in TableEl For a good fit this 
parameter should be close to unity. The likely reason that 
we have relatively high values of x^ for about half of the 
sources is that we do not take into account the uncertain¬ 
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Fig. 8. The templates used in the fitting procedure (solid 
curves) together with the resulting best fits using all other 
templates (dotted curves). The spectra are all normalized 
such that the maximum value of the template equals unity. 
In the left column we show the templates for the small 
grains (O.I/im), in the right column those of the large 
grains (1.5/rm). 


ties on the Ki. These are mainly caused by uncertainties 
in the shape, structure and size of the grains, and in the 
laboratory measurements of the wavelength dependent re¬ 
fractive indices. 

To test whether we have degeneracies between the var¬ 
ious templates used in the fitting procedure we tried to fit 
each of the silicate templates using a linear combination of 
all other templates. The results are shown in Fig. 0 where 
we plot the mass absorption coefficients together with the 
best fit using the other templates (here the mass absorp¬ 
tion coefficients have been normalized such that the max¬ 
imum value equals unity). The figure shows that almost 
all of the templates used have a unique spectral structure 
that cannot be reproduced by the other templates. Only 
the emission from large pyroxene grains can be reproduced 
reasonably well by the other templates, although signifi¬ 
cant differences still exist (e.g. the 11.3/im feature that is 
present in the fit to the large pyroxene opacity in Fig. 
The fit consists of 52% large olivine grains, 43% small 
pyroxene grains and only 5% of crystalline silicates. This 
could result in a slight change in the mass fraction of large 
grains when this template is not used. The crystalline frac¬ 
tion would not be affected significantly. Because the small 
pyroxene grains are needed in order to reproduce the short 
wavelength side of some of the spectra, we chose to include 
also the large pyroxene grains for consistency. The pres¬ 
ence of large pyroxene grains can be firmly established 
only in the highest SNR spectra. 
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Table 6. The best fit values of the parameters in our compositional fits. The abundances of small (0.1 /rm) and large (1.5 /xm) grains of the various dust species 
are given as fractions of the total dust mass, excluding the dust responsible for the continuum emission. If a species was not found, or unconstrained by the 
spectra, this is indicated by a - symbol. The PAH and continuum flux contributions (the last two columns) are listed as percentages of the total integrated flux 
over the 10/xm region, contained in these components. These are measures for the relative flux contributions, but cannot be interpreted as relative dust masses. 


# 

Star 


Tc[K] 

Olivine 

Pyroxene 

Forsterite 

Enstatite 

Silica 

PAH 

Cont. 





Small 

Large 

Small 

Large 

Small 

Large 

Small 

Large 

Small 

Large 

contr. 

contr. 

0 

Gal. Center 

66.8 

- 

85.8l?:l 

- 

12.51?;? 

- 

o.6l?;? 

- 

- 

- 

- 

4-‘4-0.0 

'^•^- 0.0 

- 

1 

AB Aur 

2.5 

464t^ 

49.511? 

47.2111 

- 

- 

0.6l?;? 

r \ y+0.7 
-0.7 

- 

- 

- 

- 

1 Q+ 0-1 
4-y_o.2 

53. il?;? 

2 

UX Ori 

0.7 


63.61?;? 

30.71?;? 

- 

2 .OI?;? 

-1 q+ 0.4 
l-'^-0.4 

o.3l?;? 

o.6l?;? 

1 1+1'^ 

- 


0 1+0-3 

45 . 21 ?;? 

3 

HD 36112 

1.7 

6691^? 

50.51?;? 

22.41?;? 

- 

22 . 2 I?;? 

2.9l?;? 

o.ol?;? 

i.8l?;? 

'-'•^- 0.2 

- 

- 

n i+'^-^ 

33.91?;? 

4 

HK Ori 

0.4 


7.11?;? 

0.9lJ®9' 

13.71?;? 

7^ q+9.9 

t 0 . 0 _i 3 3 

n 7 +^ *3 

1 4+2'^ 

4-+-1.3 

o.il?:? 

0 1 

u.r-o.i 

- 

o.stj'i 

9 p:+0.4 
^•^-0.4 

68.51?;? 

5 

HD 245185 

0.6 


28.51?;? 

57.41?;? 

o.il?:? 

12 . 2 I?;? 

0 2^^'^ 
'^“^-0.2 

Q 7 + 0.9 
‘4''-0.6 

O.Ol?;? 

0 0 +°'^ 
'-’•^- 0.0 

- 

0 0+0-5 
^•^- 0.5 

1 1 + 0.2 
4*4_o.3 

50.61?;? 

6 

V380 Ori 

1.8 

48lti 

- 

0 . 8 l?;? 

- 

74.61?;? 

Q Q + 0.6 

^•^-0.5 

10 . 8 l?;i 

- 

10 . 81 ?;? 

0 . 8 l?;? 

n 0+’^'^ 

n 7+0-1 
^ -0.1 

84.91?;? 

7 

HD 37357 

0.5 

5081"" 

0 0+5.9 
'^•''-3.4 

39.41“;? 

O.ll?;? 

46.0l?V 

Q y+0.5 
^•'-0.6 

Q q+ 1.8 

^•^-1.5 

0.9l?;? 

1.3l?;? 

O.Ol?;? 

q r+ 0.8 
'^•^- 0.8 

0 1 +'^'^ 
O'-4-0.1 

47.51?;? 

8 

HD 37806 

1.8 

5251^ 

- 

59.81?;? 

- 

0 - 8 l?;? 


3.51?;? 

O. 2 I?;? 

20.91?;? 

5.01?;? 

q r+ 0.8 
'^•^- 0.8 

i. 8 l?:? 

7q 7+1.0 

i 0. t _Q 9 

9 

HD 95881 

3.7 

42919 

- 

- 

- 

79.71?;? 

Q 7+0.5 

4.7l?;? 

0 n+°.® 

U.U-O.O 

ii.3l?:I 

3 . 0 I?;? 

n ^+00 

^•^-0.5 

3 . 0 I?;? 

84.21?;? 

10 

HD 97048 

1.4 

3031^ 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

14. ll?;? 

83. ll?;? 

11 

HD 100453 

2.1 

429 I 9 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

4.4l?;? 

93.51?;? 

12 

HD 100546 

5.6 

2611? 

2.5l?;? 

78.51?;? 

- 

o.ol?;? 

fi 2+^-^ 

o.il?:? 

- 

c 7+0.5 

' -0.4 

0 . 2 I?;? 

a 7+0-3 

q 7+0.1 
• -0.7 

55.31?;? 

13 

HD 101412 

0.5 

4251? 

- 

0 

'^•'^-0.3 

0 . 2 I?;? 

75.21?;? 

- 0.6 

3.4l?;? 

0 1 

16.81?;? 

q Q + 0.6 

o.^_0.6 

0 1 +^'^ 

Z.U-o 2 

80.21?;? 

14 

HD 104237 

3.9 

4461? 

- 

- 

- 

77.71?;? 

5 0+0-4 
o.y_o.6 

- 

- 

11.81?;? 

1 7+0.2 
4-'-0.1 

2 0+0-4 
^•^-0.4 

1 4+0.2 

4-^-0.1 

63.41?;? 

15 

HD 135344 

0.8 

680l(? 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

4 0+0-3 
^•^-0.3 

90.31?;? 

16 

HD 139614 

0.6 

2771? 

I 6 . 0 I?;? 


- 

OQ 7+13-8 

2 4+0.9 
^•^-0.8 

0 7+^ ’^ 
‘4''-0.6 

0 0 +°'^ 

'-’•V-O.O 

4 . 11 ^ 2 :? 

- 

n 1+°'^ 

9 9+0.3 
^•^-0.3 

75.91?;? 

17 

HD 142666 

1.6 

3961? 

42.81?;? 

Q 7+5-0 

- 

43.51?;? 

i-+?:3 

r \ q + 0.8 

^•O_0.8 

O.Ol?;? 

n fi+o-8 

- 

- 

1 -^+0.2 

4-'-0.1 

66.81?;? 

18 

HD 142527 

11.1 

5231? 

- 

'^•'^-0.0 

- 

73 . 2 I?;? 

8.6l?;? 

- 

0.6l?;? 

14 2+0.8 

^•^-0.1 

i.ol?:? 

3 . 1 I?:? 

58.81?;? 

19 

HD 144432 

3.0 

40ll?° 

52.71?;? 

o-ol?;? 

- 

42.31?;? 

1 Q+0.1 

Q 7+0.3 
‘4'^-0.3 

0 7+0-2 

O’*' -0.2 

n 0+0-3 

y-y-0.3 

- 

o.6l?;? 

0.51?;? 

39.71?;? 

20 

HD 144668 

9.6 

5171? 

- 

56.9l?3«2 

- 

12.0l?l^ 

7 n+0.3 

^•11-0.3 

O+0.9 

'3''^-0.8 

O.Ol?;? 

1 2 1 +'^'® 

1 k+0.2 

4*^-0.3 

4 . 2 I?;? 

1 2+°-i 

4-^-0.1 

85.31?;? 

21 

HD 150193 

5.1 

410l?o 

- 

0 

^•'-0.6 

- 

82.ll?;? 

q q+0.1 
'^•'^-0.1 

4 1 +0-3 

0 . 2 I?;? 

q 7+0.3 
' -0.4 

- 

5.81?;? 

0 2+*^'^ 
44*^-0.1 

39.71?;? 

22 

HD 163296 

10.7 

46ll? 

15.3l?;I 

29.81?;? 

- 

42.21?;? 

0 1+0.1 
'J'l-O.l 

1 q+0.3 
i--'^-0.3 

0 . 2 I?;? 

4 0+0-3 
^•11-0.3 

- 

4 . 1 I?:? 


49.61?;? 

23 

HD 169142 

0.9 

3481? 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

16. ll?;? 

82.01?;? 

24 

HD 179218 

11.8 

3501?4 

- 

1 1 Q+2-6 

^^•'^-2.6 

- 


- 

q r+0.6 
O . 0 _Q Y 

6.2l?;? 

18.31?;? 

2 n+oi 

- 

4.4l?;? 

64.81?;? 
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There are a few points one has to keep in mind when 
interpreting the results of the analysis presented above. 

— Since we consider only the 10 micron region of the 
infrared spectrum our data are sensitive only to grains 
with a temperature of ^ 200 K or more. This limits 
our study to the inner disk regions (i.e. < 10-20 AU). 

— Because of the same limitation of the spectral range 
our data are only sensitive to relatively modest sized 
grains. In order for the emission to show an observable 
spectral signature, the dust grains have to be small 
compared to the wavelength of radiation. 10 /im mea¬ 
surements are sensitive to grains with a volume equiv¬ 
alent radius rv ^m. 

— The disk regions that we study are partially optically 
thick at 10 /im. This means that we cannot see deep 
into the disk. We only observe the surface layer. 

— Due to the limited spatial resolution of our observa¬ 
tions, we observe the integrated spectra of the entire 
inner disk surface. The observed flux is therefore an 
average over distance to the star and thus over tem¬ 
perature. Also, we know from spatially resolved obser¬ 
vations of the innermost part of a few disks that the 
mineralo gy is not constant as a fu nction of distance to 
the star I van Boekel et al.l2004afl . Close to the star we 
have higher temperatures and densities which trigger 
both crystallization and grain growth. The derived pa¬ 
rameters therefore represent an average, characteristic 
temperature and an average dust composition. 


5.3. Observed trends in the fits 

We will now discuss the trends and correlations observed 
in the derived fit parameters. 

In Fig. 13 we visualize the grain growth and crys¬ 
tallinity, as implied by our compositional hts. Horizontally 
we plot the mass in large (1.5/tm) grains, as a fraction 
of the total dust mass, excluding the dust responsible for 
the continuum component, and the PAHs. Vertically, we 
likewise plot the mass fraction contained in crystalline sil¬ 
icates in small and large grains, which is also referred to 
as the crystallinity of the material. 

Upon inspection of the figure it is clear that all disks 
show signs of substantial removal of small grains. There 
are no sources with a mass fraction in large grains below 
30%. This infers that none of the sources in our sample 
contains truly “pristine” dust. All sources have an appre¬ 
ciable amount of large grains at their disk surface, com¬ 
pared to ISM conditions. In addition, all sources have a 
crystallinity that is higher tha n the value we derive for 
the ISM (~0.6%; Kemper et al. l2004l derive an even more 
stringent upper limit of 0.4%). 

The derived mass fraction in large grains ranges from 
w30% to wl00%, with most sources at high values. All 
sources exhibiting a high degree of crystallinity have a high 
mass fraction in large grains. There are no highly crys¬ 
talline sources (crystallinity above 10%) with less than 
85% of the dust mass in large grains. 



Fig. 10. The fraction of the dust mass contained in silica 
grains vs. the mass fraction cont ained in forsterite grains. 
Following lBoiiwman et all ll200lll we also plot the theoret¬ 
ical annealing behavior of two different amorphous magne¬ 
sium silicates, smectite dehydroxy late (SMD; Mg6Si8022; 
upper dashed line) and serpentine dehydroxylate (SD; 
Mg 3 Si 207 ; lower dashed line). The solid line represents 
the expected annealing behavior of a mixture of these two 
silicates, consisting of 4% SMD and 96% SD, which was 
found by iBouwman et alJ ll200lll to give the best fit to 
their data. Our data are in better agreement with a pure 
SMD initial composition. 


There are no sources with a mass fraction of crystalline 
material above 35% (see also TableEJ- It should be kept in 
mind that the silicate emission we see likely ori ginate s in 
the surface layer of the disk. Van Boekel et al. 11200311 ar¬ 
gue that the disks are well mixed in the vertical direction, 
and that therefore the observed silicate emission should be 
representative of the whole micron and sub-micron sized 
dust population of the disk. We point out that in this 
work, “crystallization” refers to the process of crystalliz¬ 
ing the material (by whatever means), spreading it over 
a significant part of the disk region seen at 10 /tm, and 
bringing it up to the disk surface where we can see it 
spectroscopically. In a scenario where the crystalline sili¬ 
cates are produced by thermal annealing in the innermost 
disk regions, and transported outward by radial mixing, 
the degree of crystallization is therefore a measure of the 
degree of mixing in the disk rather than the actual pro¬ 
cess of annealing, which is effectively instantaneous at the 
inner disk edge. 

In Fig. cni we plot the mass fraction of dust con¬ 
tained in forsterite versus that contained in silica grains. 
It is clear that these mass fractions are correlated. 
Experiments show that when forsterite is created by an¬ 
nealing o f an amorphous silicate, silica is formed as a by¬ 
product llRietmeiier et al.lll986tlHallenbeck fc Nuthlll997t 
iFabiau et a.lJ 1200611 . The amount of silica that is created 
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Fig. 9. The relation between grain size and crystallinity found in our spectral fits. Vertically, we plot the mass fraction 
of crystalline grains (forsterite and enstatite). Horizontally, the mass fraction of large (1.5/im) grains is plotted. Since 
groupIb sources do not display a silicate feature, they are not shown here. For a discussion see sections l5.dl and El 


when forming a certain amount of forsterite depends on 
the type of amorph ous silicate one starts o ut with. As al¬ 
ready suggested bv iBouwman et all l 200lll . we can try to 
constrain the composition of the amorphous material by 
measuring the ratio of forsterite over silica. In the figure 
we also plot the expected annealing behaviour of smectite 
dehydroxylate (SMD; Mg6Si8 022) and serpentine dehy- 
droxylate (SD; Mg 3 Si 2 07 ). In low temperature conden¬ 
sation experiments, th ese are the only magn esium sili- 
cates that are f o rmed llHietmeiier et al Whereas 

iBouwman et al ] (l 2 nmh found that a mixture of 4% of 
SMD and 96% of SD yielded the best fit to their data, 
our results seem to favour an initial composition of pure 
S MD. Possibly, the discrep ancy between the results found 
by iRoiiwman et all ll2nfl1ri and the results found here is 
connected to the differences in th e dust componen t s used 
to fit the 10/xm spectra. Notably. iBouwman et al ] ll20nih 
do not include large silica and forsterite particles. 


In Fig. we show the mass fraction of crystalline ma¬ 
terial contained in enstatite versus the total mass fraction 
of crystalline material. There is a correlation between the 
fraction of the total mass contained in crystalline silicates 
and the composition of these crystalline silicates. In gen¬ 
eral, for sources with a high degree of crystallinity most 
crystals are in the form of enstatite, while for the sources 
with a low crystallinity, forsterite is the dominant crys¬ 
talline species. We will discuss this further in section ESI 

The amount of growth that the crystalline material has 
experienced is compared to the growth in the amorphous 
component in Fig. ESI If the amorphous grains are large, 
also the crystalline grains are large, though the correlation 
is not tight. In all sources in which the amorphous compo¬ 
nent has more than 85% large grains, also the crystalline 
component is dominated by large grains. The sources that 
have less than 85% large grains in the amorphous com¬ 
ponent all have a low crystallinity (see Fig. El). Therefore, 
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Fig. 11. The mass fraction of crystalline material con¬ 
tained in enstatite vs. the total mass fraction of crystalline 
material. A value of 0 on the vertical axis indicates that 
all the crystalline silicates present in the disk are in the 
form of forsterite, while a value of 1 means all crystalline 
silicates are in the form of enstatite. 



Mass Froction in Large Grains in Amarphous Component 


Fig. 12. The mass fraction of large grains in the crys¬ 
talline grain population vs. the mass fraction of large 
grains in the amorphous component. When the average 
size of the amorphous grains is large, the majority of the 
crystalline material resides in large grains as well. 


the ratio of large and small crystals is poorly constrained 
in these sources, which is reflected in the large errorbars. 

The fraction of crystalline silicates is correlated with 
the mass and luminosity of the central star. This is vi¬ 
sualized in Fig. El The higher mass (higher luminosity) 
stars have an on average higher crystallinity than the lower 



Fig. 13. The mass fraction of dust in crystalline grains vs. 
the stellar mass. Higher mass stars show an on average 
higher fraction of crystalline grains than do lower mass 
stars. 



Fig. 14. The mass fraction of dust residing in large grains 
vs. the stellar mass. All stars with a mass M > 2.5 M© 
have a mass fraction of large grains above 85 %. 

mass (lower luminosity) stars. This correlation will be dis¬ 
cussed in section lO 

All stars with a stellar mass above about 2.5 M© 
{L >60Lq) have a high fraction (>85%) of large grains 
(Fig. [nil . Possibly, conditions in the disks around more 
massive stars are more favourable for growth than they are 
in the disks around lower mass stars. An alternative expla¬ 
nation for the observed trend is that the disks around the 
more massive stars, which are all relatively young (Fig.|21), 
are more turbulent than those around the less massive, 
older stars. As a consequence the mixing in the young 
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disks will be more efficient. While in the older stars the 
large grains decouple from the gas and settle to the mid¬ 
plane (and therefore will not be detected in the 10 micron 
spectrum anymore), in the young stars, larger grains may 
still reach the disk surface and cause the observed average 
grain size to be higher. 

6. Discussion 

We now discuss our fit results in terms of the processes 
that are responsible for the dust evolution. We will first 
briefly outline the expected conditions that prevailed dur¬ 
ing the active disk phase ('section liTTll . i.e. the phase prior 
to the passive disk phase. For a compre hensive review of 
active disks, see e.g. ICalvet et all ll2n0(Tll . The disks of all 
stars in our sample are in the passive disk phase. The 
characteristics of this phase are discussed in section liT^ 
In sectio nifi.dl we summarize the constraints put on dust 
processing by this work and previous studies. Lastly, in 
sectio ulb.di we sketch a scenario that is consistent with the 
current knowledge of dust processing and disk evolution. 

6.1. The active disk phase 

The accretion of matter onto a forming proto-star is be¬ 
lieved to be initially spherically symmetric. At some point, 
as the proto-star contracts, conservation of angular mo¬ 
mentum inhibits further spherical accretion. The accretion 
process then proceeds through a disk, and is accompanied 
by a bi-polar outflow. The fact that material is radially 
transported through the disk implies that the disk has 
viscosity, which is in turn coupled to turbulence. 

During this so-called active disk phase, gravitational 
energy of the accreting material is dissipated in the disk, 
thereby heating it. Accretion rates 
may reach values up to 10 “^ M 0 yr“^ 

Close to the central star the main energy source of the 
disk is accretion luminosity, whereas at larger radii irra¬ 
diation by the star and hot inner disk regions is expected 
to be the main heating source of the disk material. 

The dust in the disk consists mainly of silicates (stud¬ 
ied in this work), and carbon. The dust does not con¬ 
tribute significantly to the disk mass in this phase of the 
disk evolution (the gas over dust ratio is on the order of 
10^ by mass). However, the thermal radiation emitted by 
the dust is the dominant cooling process in the disk. In 
the outer disk regions where the main energy source is 
irradiation, the dust governs the heating of the disk as 
well. Therefore, the dust properties determine the disk 
temperature, except in the innermost disk region where 
the heating is dominated by viscous dissipation of gravi¬ 
tational energy of the accreting material. The gas, heated 
by the dust, provides the pressure support of the disk. 

The densities in the disk are high, especially near the 
disk mid-plane, and one may expect coagulation of small 
dust grains into larger aggregates to occur. Already after 
~ 10 ^ yr of dynamical disk evolution the average mass of a 
grain can increase by a factor 10 ^ to 10 ^, and close to the 


in the active phase 



disk midplane the grains may re ach sizes of at least a few 
times 10/im up to millimeters llSiittner fc Yorkj l200ltl . 
Close to the central star the temperatures in the disk 
can reach values in excess of 1 000 K. At this temperature 
the (initially mostly amorphous) silicates are annealed, to 
form crystalline silicates. Yet closer to the star, temper¬ 
atures reach values of about 1500K, and crystalline sili¬ 
cates may form by gas phase condensation of evaporated 
material. 


In some of the stars in our sample, the crystalline sil¬ 
icates appear so prominent that it is unlikely that their 
emission arises from the innermost disk region only; the 
abundance of crystalline silicates must be high in a rel¬ 
atively large part of the region of the disk responsible 
for the 10 ^m emission (the innermost 10-20 AU of the 
disk). There are in essence two possible ways to get crys¬ 
talline material at ~10AU distance from the star: (1) 
thermal processing in the hot inner disk and subsequent 
radial transport of this material outward, and ( 2 ) local 
production of crystalline material at large distance from 
the star in transient h eating events cau sed by shocks (e. 
Harkei_^ Desc]J 20 (i^ or lightning (e.g. IPilioD et al.ll99 


Desch fc Cuzzi200(]ll. During the active disk phase the ac¬ 


cretion rate in the viscous disk is high and the disk will 
be turbul ent. This is expected to enable radial m ixing of 
material I Bockelee-Morvan et al 1 l2002HGaillE00l . 


6.2. The passive disk phase 


When the supply of fresh material from the maternal cloud 
has exhausted, further accretion onto the star proceeds 
only on a very low level (order 10“® M 0 yr“^). Under these 
conditions the energy production by viscous dissipation 
can be fully neglected. Throughout the disk the tempera¬ 
ture of the material is determined by absorption of stellar 
radiation; this is referred to as the passive disk phase. The 
stars that we study are in this evolutionary phase. In the 
observed SEDs of sources in the passive disk phase, the in¬ 
frared excess typically becomes noticeable above the pho- 
tospheric emission at a wavelength of about 1 /im. The ex¬ 
cess emission in the near-infrared indicates that the inner 
radius of the disk is determined by the evaporation tem¬ 
perature of silicates (at about 1500K). At wavelengths 
longward of 2 /im the infrared emission from a gas rich 
disk completely dominates the SED. 

The formation of planets and planetesimals is thought 
to occur during the passive disk phase. As the disk dissi¬ 
pates on a timescale of 10 ^ yr, the infrared excess fades. 
The inner disk regions become devoid of gas and dust first. 
When the syst em has evolved into the debris disk p hase 
(e.g. 3 Pictoris. lAuman iil llflSdt ISmith fc TerrilelllQS^ ex¬ 
cess emission can be seen above the stellar photospheric 
emission only at wavelengths above ~I 0 /tm. 
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6.2.1. The era of crystallization 


During the passive disk phase there is a region in the inner 
disk where the temperature is above 1 000 K, and therefore 
the silicates will be crystallized. As there is no significant 
accretion luminosity, this region will be smaller than in the 
active disk phase. It is unclear whether the disk will be tur¬ 
bulent enough for significant radial mixing to take place 
during the passive phase. Also the proposed mechanisms 
for local production of crystalline silicates in the outer 
disk regions are more efficient in the active disk phase. 
Theoretically therefore, the active disk phase is the pre¬ 
ferred era for the crystallization of the dust. Nonetheless 
it has, from an observational point of view, been suggested 
that the crystallinity of HAe disks gradually evolves from 
low t o high values during the passive phase I Gradv et al.l 

l2nn(t . 


6.3. New constraints on dust processing 

We first briefly repeat the results that follow directly from 

our spectral modeling: 

1 . All disks have already substantial removal of the small¬ 
est grains (in our fits represented by the 0.1 pm grains). 
The mass fraction of large (in our fits 1.5/rm) grains 
ranges from about 35 to almost 100 percent. 

2. The crystallinity ranges from about 5 to 30 percent. 
Disks with a high crystallinity are always dominated 
by large grains. 

3. Vice versa, disks with more than 85 percent large 
grains show a crystallinity above 10 percent; disks with 
less than 85 percent large grains show a crystallinity 
below 10 percent. 

4. Large crystals are found when large amorphous grains 
are abundant. We do not find disks in which the bulk 
of the grains is large, and the crystalline silicates are 
small. 

5. In disks with a crystallinity below 15 percent, more 
than 45 percent of the crystalline silicates consist of 
forsterite. Above 15 percent crystallinity, more than 45 
percent of the crystalline silicates consist of enstatite. 

When the stellar parameters of our sample stars (see 

Tabled are taken into account, several additional con¬ 
clusions can be drawn: 

6 . The disks with the highest degree of crystalline ma¬ 
terial (between 20 and 30 percent) all belong to stars 
with a mass above 2.5 Mq, and a luminosity above 
60 Lq. 

7. Below a crystallinity of 20 percent, there does not 
appear to be a correlation between stellar parame¬ 
ters (mass, luminosity, age) and crystallinity. Separate 
studies of the 10 micro n spectra of TTauri stars 
(Przygodda et al. 120031: Honda, private communi¬ 
cation) indicate that the crystallinities of T Tauri 
star disks are similar to those observed in the 
disks of those HAe stars in our sample that have 


M <2.5M0 (i.e. <10 percent). We note that a small 
subset of the TTauri star disks shows a significantly 
hig^er^crystalfinitv, deviating from the general trend 

8 . The stars in our sample with masses below 2.5 Mq have 
ages ranging from 2x10® to ~10^ yr (see Fig. 0. If 
we consider only this subgroup, we find no relation 
between the stellar age and the disk crystallinity. 

9. The stars that lack a 10 micron silicate emission band 
(HD 97048, HD 100453, HD 135344 and HD 169142) 
tend to be old, with ages of 5-10 Myr. 

A number of previous studies have been concerned with 
the crystalline silicates in HAe star disks. Here we briefly 
recall some of the results that are of relevance for this 
study: 


10. Recent observations usi ng the 10/rm instrument 
MIDI llLeinert et alJ 200^ on the VLT Interferometer 
dGIindemaniT'et^ J|2Q0^ have enabled the extraction 
of the spectra of the innermo st ~2AU of three HAe 
stars l|va.n Boekel et '^ l2004all . showing that: 

- HD 144432 has virtually all its crystals in the inner 
2AU. 


- HD 163296 has a flatter crystallinity ’’gradient” and 
has a higher overall crystallinity than HD 144432. 

- HD 142527 has the highest crystallinity of the three 
stars studied, virtually the entire inner disk is crys¬ 
tallized. Also at radii larger than 2 AU, the abun¬ 
dance of crystalline silicates is relatively high. 

- The innermost regions of the disks have experi¬ 
enced more grain growth than the outer disk re¬ 
gions. Small grains are mostly found at larger dis¬ 
tance from these stars. 

11. ISO observations of HD 179218 and HD 100546 show 
the presence of cold (100-150 K) crystalline silicates. 
This implies a substantial cryst allinity even at dis¬ 
tances between 20 and 40-50 AU llMalfait et alJll998b[ 
iBouwman et al.ll200lL MEOl). 

12. A 33.5 pm forsterite band is te ntatively detected in 
the IS O spectrum of HD 100453 llVandenbussche et alJ 
I2flfl4ll . 


We will now discuss the implications of the points men¬ 
tioned above. For clarity, we have labeled each of the im¬ 
plications with the index number(s) of the point(s) that 
lead to it. 

1, 2, 3: The dust in the circumstellar disks coagulates 
more easily than it crystallizes. This is likely due to 
the fact that the circumstances enabling grain growth 
(i.e. high densities) prevail in a much larger part of 
the disk than the circumstances needed for crystalliza¬ 
tion of silicates (i.e. high temperatures). Once a certain 
level of coagulation is reached, crystallinity becomes a 
lot more obvious in the upper layers of the disk. 

4: This observation is consistent with crystallization oc¬ 
curring as the grains grow: whatever population of 
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grains is present, is crystallized. Most of the crystal¬ 
lization occurs when the average grain size is already 
large (micron sized). 

1, 2, 3, 10, 12: Coagulation is more efficient in the inner 
disk regions than in the outer disk regions. In the 
sources that lack a silicate band coagulation must 
have proceeded furthest. Our relation between grain 
growth and crystallization (Fig. then implies that 
these disks must have a substantial fraction of crys¬ 
talline silicates, certainly above 10 percent and per¬ 
haps more than 20 percent. However, we do not find 
the usual forsterite band at 11.3/rm nor the 23.5/rm 
forsterite band. So, also the crystalline silicates must 
on average be l arge. Indeed this is what was concluded 
for HD 100453 llVandenbussche et a i] |2004ll . Again, it 
seems that whatever (inner disk) grain population is 
present, is crystallized. 

4: Micron sized crystalline grains are present in the disks. 
It is hard to produce crystalline grains of these sizes 
via shock heating. At any rate, our data put severe 
limits on the efficiency of shock heating, as the shock 
mechanism must be able to produce sufficiently strong 
shocks to crystallize micron sized grains. 

5: If we consider two possible sources of crystalline sil¬ 
icates, i.e. chemical equilibrium reactions (gas phase 
condensation and subsequent gas-solid reactions) and 
thermal an nealing, then we can note the following 
llGa,illl2nnl : 


- Enstatite is expected to be the dominant crystalline 
species formed by chemical equilibrium processes 
in most of the inner disk, apart from the hottest, 
innermost region, which is dominated by forsterite. 

- Since most likely amorphous olivine is the most 
abundant species entering the disk from the ISM, 
forsterite (or better, crystalline olivine) is expected 
to be the dominant crystalline species formed by 
thermal annealing. 

In disks with a low crystallinity, forsterite is the dom¬ 
inant crystalline species, suggesting that this material 
is formed by thermal annealing. The high abundance 
of enstatite observed in the highly crystalline disks in¬ 
dicates that the production of crystalline material has 
occurred by means of chemical equilibrium processes 
in these sources. 

I, 2, 3, 4, 10: The TIMMI2 and MIDI data taken together 
are consistent with crystallization starting in the in¬ 
nermost regions. Disks with a high crystallinity in the 
TIMMI2 data have a substantial fraction of crystalline 
grains in the 5-10 AU area (or even in the 20-50 AU re- 
gion. see e.g. the ISO spectra of HD 100546, H D 179218 
(jMalfait et fil.l[T998bl:lRoiiwmfi,n et fil.ll200l|B . i.e. well 
outside the regions where thermal annealing and chem¬ 
ical equilibrium processes are expected to be effective. 

5, 11: HD 179218 has a high abundance of cold enstatite. 
In the context of the above reasoning, this enstatite 
must have been produced either locally in transient 
heating events, or transported outward from the inner 
disk to observable cold disk regions. If chemical equi¬ 


librium processes are an important means of producing 
enstatite, the HD 179218 data point to enstatite pro¬ 
duction in the inner regions and transport to larger 
distance by radial mixing. 

6: Crystallization is more efficient in disks around higher 
mass, higher luminosity stars. Several possible expla¬ 
nations for this trend can be identified: 

1. during the active phase, the region in which the 
disk material is sufficiently hot for thermal anneal¬ 
ing (and chemical equilibrium processes) to occur 
is larger in high mass stars, since the energy dissi¬ 
pated in the disk per unit mass accreted material 
is higher. If this is the dominant effect, our results 
favour crystallization in the active phase. 

2. during the passive phase, when the extent of the 
region where annealing (and chemical equilibrium 
processes) can occur is set by irradiation from the 
central star, the higher luminosity of the more mas¬ 
sive stars causes this region to be larger than in the 
low mass stars. If this is the dominant effect, our 
results favour crystallization in the passive phase. 

3. the disks around more massive, more luminous 
stars may be more turbulent, both during the ac¬ 
tive and passive phase. The outward mixing of pro¬ 
cessed material will then be more efficient in the 
high mass stars, causing a more prominent appear¬ 
ance of this dust in the spectrum. 

It is currently unclear which of the above effects is 
dominant. The study of young (r < 10® yr), low mass 
stars may allow to answer this question. 

8: This observation suggests that whatever caused the ob¬ 
served range of crystallinities in stars with a mass be¬ 
low about 2.5 Mq, occurred before 2x 10® yr. This puts 
the epoch of crystallization in the active and/or early 
passive phase. 


6.4. Global picture that emerges 

The 10 pm emission is emitted by the surface layer of the 
disk at radii below 20 AU from the central star. In the mid¬ 
plane of this region (and presumably even at larger radii), 
dust coagulation takes place. Vertical mixing ensures that 
the coagulated, micron-sized grains reach the disk surface 
layer and become apparent in the spectrum. Growth is 
most efficient in the innermost disk regions where the den¬ 
sities are highest. Crystallization by means of thermal an¬ 
nealing and chemical equilibrium processes occurs in the 
hot innermost disk region. In many disks there are signifi¬ 
cant amounts of crystalline material at larger (>5 AU) dis¬ 
tances from the star, where thermal annealing and chem¬ 
ical equilibrium processes are ineffective. The combined 
evidence presented in section lO seems to favour radial 
transport from the innermost disk regions as the source 
of these crystals, above local production mechanisms (e.g. 
shock annealing). 
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Crystallization is most efficient in the disks surround¬ 
ing the more massive, more luminous stars (M> 2.5 Mq, 
L> 6 OL 0 ). In these disks, the region in which thermal 
annealing and chemical equilibrium processes can pro¬ 
duce crystalline silicates is larger than in the disks sur¬ 
rounding lower mass stars. Additionally, the disks around 
the more massive stars may be more turbulent, enabling 
more efficient radial mixing. It is unclear whether the 
crystallization of the dust in the disks occurs predomi¬ 
nantly in the active disk phase, or in the passive phase 
that follows. Our data do suggest that the disks reach 
their final crystallinity relatively early in the passive phase 
(t <2 Myr). Therefore, crystallization happens during the 
active and/or early passive phase. 


7. Conclusions 


We have undertaken a large spectroscopic survey of Herbig 
Ae stars in the 10 /rm atmospheric window, and have pre¬ 
sented new spectra of 23 stars. The correlation between 
the shape and the str ength of the 10 um silic ate emis¬ 
sion band reported bv Ivan Boekel et ^ 1 2003h is recon¬ 
firmed with a larger sample. We have performed composi¬ 
tional fits to the silicate feature using opacities of minerals 
commonly found in circumstellar disks (amorphous olivine 
and pyroxene, crystalline olivine and pyroxene, and amor¬ 
phous silica) and polycyclic aromatic hydrocarbons. For 
all minerals we have allowed both “small” ( 0.1 ^m) and 
“large” (1.5 /im) grains in order to study the effects of 
grain growth. This set of materials is sufficient to repro¬ 
duce the wide variety of observed spectral shapes. 

We find a trend between the mass fraction in large 
grains and the mass fraction in crystalline grains: all 
sources with a high crystallinity have a high mass fraction 
in large grains. There are no highly crystalline sources 
which are dominated by small grains. Most sources have 
a mass fraction in large grains of more than 80 percent. 

We note that there is an important bias in our sam¬ 
ple (and probably in most samples of Herbig stars studied 
in the literature): the more massive sample stars (~ 3 M 0 , 
“high mass”) are younger than the less massive (< 2.5 M 0 , 
“low mass”) stars. There are no low mass stars younger 
than 1 Myr in our sample. In order to establish the disk 
conditions at the end of the active disk phase, which pre¬ 
cedes the passive phase, it is essential that such very 
young, low mass stars are found and studied. Since these 
stars may still be enshrouded in circumstellar material, 
they may have to be selected using infrared data. 

We find a trend between the derived crystallinity of the 
dust and the mass (and luminosity) of the central star: the 
disks around stars with a mass larger than 2.5 M 0 (a lu¬ 
minosity above 60 L 0 ) have a higher crystallinity (> 20 %) 
than the less massive, less luminous stars. Within the sub¬ 
set of sources with a stellar mass below 2.5 M 0 , no corre¬ 
lation between crystallinity and stellar parameters (mass, 
luminosity, age) is seen. These lower mass stars in our 
sample are all older than ~ 2 Myr. Since in this subset 
there is no correlation between age and crystallinity, we 


conclude that the crystallization of the material predom¬ 
inantly happens in the active or early passive disk phase 
(before 2 Myr). 

The evidence presented in this paper combined with 
conclusions from other studies seems to favour a scenario 
in which crystalline silicates are produced in the innermost 
regions of the disk and transported outwards. Spatially 
resolved spectra of these disks, as can be obtained us¬ 
ing for example the MIDI instrument on the Very Large 
Telescope Interferometer, will provide crucial information 
on the radial dependence of the mineralogy of the dust 
in these disks. In addition, measurements at longer wave¬ 
lengths can probe colder regions (further out) in the disk 
which can provide further constraints on the temperature 
structure and spatial distribution of the dust. 
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